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Abstract 

The group velocity of gravitational waves in a flat Friedman-Robertson- Walker universe is 
investigated. For plane waves with wavelength well inside the horizon, and a universe fllled with 
an ideal fluid with the pressure to density ratio less than 1/3, the group velocity is greater than 
the velocity of light. As a result, a planar pulse of gravitational waves propagating through the 
universe during the matter/dark energy dominated era arrives to the observer with the peak 
shifted towards the forefront. For gravitational waves emitted by inspiralling supermassive black 
holes at the edge of the observable universe, the typical shift that remains after the effects of 
nonplanarity are suppressed is of order of ten picoseconds. 

1 Introduction 

Recent efforts in the construction of large detectors of gravitational waves based on laser inter- 
ferometry (for a review, see [ij) revived interest in investigations of emission, propagation and 
detection of gravitational waves. Topics that have been studied in some detail in the past in- 
clude propagation of weak gravitational waves in a Friedman- Robertson-Walker (FRW) universe 
[nisi mis]. The main result is that the waves do not obey Huygens' principle unless the universe 
is fllled with pure radiation (so that its scalar curvature is zero). Otherwise there appears a 'tail' 
in the Green's function, coming from the backscattering of waves on the spacetime curvature. 
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This contrasts with the behavior of the electromagnetic waves, whose 'tail', if present, consists 
of pure gauge. 

A kinematic characteristic often used in the description of wave propagation is the group 
velocity. Strictly speaking, it is defined for modulated waves only, as the velocity with which 
'wave groups' travel through space. However, one commonly speaks of the group velocity of wave 
pulses of limited duration, too, referring to the velocity of the peak of the envelope of the pulse. 
(For compact pulses, this can be identified with the signal velocity, defined as the velocity with 
which the main part of the signal arrives to the observer [6].) Here we address the question what 
is the group velocity of gravitational waves in a flat FRW universe. We restrict ourselves to a 
flat 3-geometry because it is easier to deal with than a closed or open one, and is preferred by 
observations. From the results on the validity of Huygens' principle it follows that if the universe 
is filled with pure radiation, plane gravitational waves advance as a whole with the velocity of 
light (along null geodesies with respect to the background metric) . In this case the group velocity 
is, of course, also the velocity of light. For other kinds of cosmic environment the group velocity 
is presumably different; moreover, while the 'tail' of the Green's function lays entirely inside the 
light cone, there is no a priori reason for the 4- vector of the group velocity to lay inside it, too. A 
superluminal group velocity does not imply violation of causality since the information carried by 
a pulse arrives to the observer with the forefront, which travels with the velocity of light (as follows 
from the Hadamard discontinuity formalism, see [7]). This is similar to what happens when light 
passes through a medium with anomalous dispersion. As demonstrated experimentally, a pulse 
of light can propagate, without being significantly distorted, even with negative group velocity 
[8]. However, no problem with causality arises since in a dispersive medium, too, the forefront of 
a pulse travels with the velocity of light [6] . 

In section [2] we compute the group velocity of gravitational waves using a well-known formula 
from optics, in section[3]we make sure that the formula can be applied to the problem under study, 
and investigate its subtleties, within the theory of wave propagation in a time-dependent medium, 
and in section 0] we discuss the results. We use the signature of metric tensor (— , +, +, +) and a 
system of units in which c = G = 1. 

2 Group velocity: calculation 

Consider weak gravitational waves propagating in a flat FRW universe. The background metric 
is 
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where t is conformal time, are comoving coordinates and a — a{t) is the scale parameter. 
Gravitational waves are described by the perturbation to the metric of the form 



a^eijdx^dx^ 

where is the traceless transversal part of the total perturbation, en — j = 0. Denote the 
derivative with respect to t by an overdote. The Hnearized Einstein equations without anisotropic 
stress yield (see, for example, [9]) 

Cij + 2Hetj - Aey = 0, 

where H ~ d/a is the Hubble parameter. After passing from to Eij — acij, one arrives at an 
equation in which the damping term is traded for a term with an external field. The equation 
reads 

- AE,j + VE,j = 0, (1) 

where 

V = -H-H^. (2) 

Suppose the universe is filled with a perfect fiuid (possibly, including dark energy) and denote 
the density of the fiuid by p and the pressure of the fiuid by p. Then, 

H^ = —p^^ H = - — {p + 3p)a\ 



so that 



V^-^iP-^P)- (3) 



For a mixture of nonrelativistic matter (p ^ p), radiation (p — p/3) and dark energy (p — —p) 
one obtains V < 0, with the equality vaHd only if the universe is filled with pure radiation. For 
Eij, we have a wave equation in flat spacetime, with a time-dependent dispersion relation 



Lu = x/k^ + V, (4) 

where uj is frequency and k is wave number (the magnitude of the wave vector) . Suppose uj is 
real, that is, there holds either V = or ^ < and A; > |y| (the wavelength is less than the 
Jeans length). Then, the field Eij is in an oscillatory regime and one can ask how fast does it 
propagate. According to the textbook formula (see, for example, |10]), the group velocity in a 
medium with an isotropic dispersion relation is w = dk(^. This yields 

k 

V = , (5) 
so that from F < we have v>l. Note that v = 1/vph, where Vph = uj/k is the phase velocity. 
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When computing v, we have ignored the fact that uj depends on t. However, the formula 
V = dki^, supposed to be vaUd for = uj{k), should stay vaHd also for lu — uj{k, t) if the period 
of the oscillations of the field is much less than the time within which uj varies considerably. 
The period is 27r/ci; and the time scale of the variation of uj is (jj/I^I, therefore ([5]) must be 
supplemented by the condition 

« 1. (6) 

After inserting here from (|4]) we find 

» \V\ + 

The universe is dominated first by radiation, then by matter, and then by dark energy, with the 
function a{t) subsequently of the form aoct, {t — to)^ and {ti —t)^^. When using these functions 
to compute V{t), we find that in the first era, as we already know, V — 0, and in the next two 
eras V ^ — At^^, where At ^ t — to and ti — t in the second and third era respectively. The 
time At approximately equals the horizon length, if 'horizon' is understood as particle horizon in 
the matter dominated era and event horizon in the dark energy dominated era. For V ~ —At^^ 
there holds V ^ ±At~^, so that the two terms in the condition for k are of the same order of 
magnitude. As a result, the condition simplifies to 

» \V\. (7) 

From V —At~^ we find k ^ At^^ (horizon length)"^; thus, in order that a wave packet 
propagates with the group velocity |[5|, its mean wavelength must lay well within the horizon. 
For such packets, the group velocity is close to 1, with the correction of the form 

A^ = -^. (8) 

Light signals (particles moving along null geodesies with respect to the background metric) 
have coordinate velocity 1. Thus, the constraint v > 1 means that the gravitational waves either 
propagate with the velocity of light or are superluminal. The former case occurs in the radia- 
tion dominated era and the latter in the matter/dark energy dominated era. When combining 
equations ^ and ^ , we find that the group velocity differs from the velocity of light by 

47r , 
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where x = k/a is the true wave number. For an estimate, suppose the universe is filled with 

incoherent dust (p = 0). Then p ~ t~^, where r is cosmological time, and the correction to v is 

of order {kt)~'^ . As a source of gravitational waves, consider a binary system consisting of two 

black holes with masses M ^ IQ^M , orbiting around each other with the period T ~ lO'^s. 

O 



(The period drops to this value about a month before the merger.) Suppose, furthermore, that 
the system is at the edge of the observable universe, at the distance / ~ 3 Gpc from Earth. 
The time it takes for a packet of gravitational waves to propagate from the binary to Earth is 
Tprop ^ T ^ 10 Gy, and the time by which the peak of the packet arrives earlier due to the 
correction to v is 

rp2 

At ~ (xT)-Vp™p - 10-"s. 

[Zn ) T 

For a universe with the actual values of cosmological parameters, Hq = 71 km s^^ Mpc^^, 
ilm — 0.26 and JIa = 0.74, and a source with the period T — 1 hour at the redshift of the most 
distant quasar, z = 6.4, a detailed calculation yields At = 6.6 x 10^^^ s. 

3 Group velocity: theory 

Consider a field u satisfying the same equation as Eij , 

il- Au + Vu = 0. (9) 

This can describe gravitational waves in an expanding universe as well as light in a homogeneous 
time-dependent optical medium. Thus, the following considerations are only loosely related to 
general relativity; instead, one can view them as a discussion of a special problem of optics. 

3.1 Waves with definite wavelength 

Choose u in the form of a plane wave with a definite wavelength propagating along the x-axis, 
u — [/(t)e*'^^. The wave number k can now assume any real value, however, for our purposes it 
is sufficient to suppose as before k > 0. From Q we obtain 

U + uj^U = 0, (10) 

with the function uj — u}{k,t) defined in |(H). If we exchange t with x and uj with k, we obtain 
an equation for monochromatic light in a static inhomogeneous medium; and if we replace, in 
addition, lu hy E (the energy of the particle) and khy p (the momentum of the particle), we obtain 
Schrodinger equation in one dimension with a modified relation between p and E, p ^ y/E'^ + V 
instead of p = ^/2{E~V). 

For light propagating in a medium with dispersion relation u = uj{k,t), the high-frequency 
condition ^ defines the geometrical optics approximation. An analogical condition is known to 
define quasiclassical, or WKB, approximation in quantum mechanics [H], and geometrical optics 
approximation in a static inhomogeneous medium in optics [12]. To obtain an approximate 
expression for C/, write U = e~^^ and expand the complex phase if} in the powers of uj~^. If we 
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denote the terms of the order minus one and zero by V'o and the terms of the first order by Ai/i, 
we have (see the corresponding formulas in [TT] ) 

/I f ' 
u^dt-Uog^, A^p^-£^--J ^^dt. (11) 

The signs in these expressions are chosen in such a way that the wave propagates in the pos- 
itive x-direction. The first expression determines the oscillating field in the geometrical optics 
approximation, and the second expression can be included into the slowly varying amplitude and 
regarded as the post-geometrical optics correction to it. 

If the dispersion relation assumes the form (JH) , the condition ((6]) for oj implies the condition ((?]) 
for k. (This holds irrespective of the form ofV{t). In general, we have also condition fc^ ^ 
however, it can be skipped if we choose V ~ —At^^.) With /c's satisfying the condition the 
dispersion relation can be written approximately as 

^ = fc+^. (12) 

The fact that we have replaced the condition for lu by the condition for k suggests that the 
concept of geometrical optics approximation must be reconsidered. The small parameter of the 
theory is {kAt)~^, where At is the time scale of the variation of V, rather than (ujSt)'^, where 
St is the time scale of the variation of w, therefore it is advisable to expand all quantities in 
the powers of fc~^ rather than In particular, the approximate expression for u) should be 

written u — loq + Alo, with the frequency in the geometrical optics approximation ujq = k and 
the post-geometrical optics correction to it Au — V/{2k). This yields 

^o = kt-^ log k, Aij^^J Vdt. (13) 

The previously obtained expression for A-ip does not contribute here since it is of order k^^. 
Note, however, that it would contribute if we used cosmological time instead of conformal, since 
then it would be of order k~^. 

3.2 Group velocity of a Gaussian packet with infinitesimal dispersion 

From plane waves with a definite wavelength one can form a planar wave packet 

oc 

,14) 

— OO 

where a is an overall amplitude peaked at the given wave number k and differing significantly 
from zero only in an interval of a size a k. Note that the integral contains modes with all 
fc's, and we have derived formulas for ij^lk, t) only for k positive and large enough to satisfy the 



6 



high-frequency condition. However, the modes with k either positive and small or negative do not 
have significant effect on the results in the limit we are interested in, therefore we can describe 
them by an extrapolated formula for large positive k. 

Define the envelope of the wave packet as the curve \u{x, t)\ at fixed t, and the group velocity 
V as the velocity with which the peak of the envelope moves along the a;-axis. That is, 

v = xm, (15) 

where Xm is the value of a; at which the function |u(a;,t)| has maximum at fixed t. To obtain an 
approximate expression for v, let us adopt two simplifications. First, consider a Gaussian packet, 

1 / AK^\ , , 

«G = ^^expl -^-2- 1 , (16) 



and second, expand the phase in the powers of AK = K — k, 

ij{K, t)=jp + ^AK + ^TjAK^ + U^K"^ + ■■■, (17) 

where = il){k,t), ^ = rj = d'^tjj, C = d^ip, ■ . ., and compute v in the linear approximation 
in which the expansion is truncated after the ^-term. 

Let U be the function obtained by factorizing out the quickly oscillating phase factor from u, 
u = e~'^^~^^''''U. For a Gaussian packet, U in the linear approximation is 

and its magnitude is 

where the indices R and / denote the real and imaginary part of the quantity in question. We 
can see that the wave packet in a;-space is Gaussian, like in A;-space, with the width ax = l/cr 
and the peak located at 

XM = ^R- (18) 

For a general time-dependent medium we have 

ioR = J dkcodt, A^R = -^dk + ^ y" '^dt 

From the first expression we obtain the leading term in v, which is just the standard group 
velocity cited in the previous section, 

vo = dkU). (19) 
The second expression yields the correction 



For a medium with the dispersion relation ([4]) we have 

which yields v — 1 + Av with Av given by equation ([8]). In this case the group velocity is given, 
within the current accuracy, entirely by the formula v — dkOJ- 

If we include higher order terms into the Aif-expansion of ijj into the calculation, we obtain 
corrections of order cr^, cr'*, ... ioU near x — xm, producing corrections of order cr^, a^, ... to 
xm- Thus, the values of xm and v obtained for a Gaussian packet in the linear approximation 
can be viewed as exact values in the limit a = (vanishing width of the packet in fc-space). 
Equivalently, one can speak of the limit ax — oo (infinite width of the packet in a;-space) . 

3.3 Effects of finite dispersion and nongaussianity 

Consider first a Gaussian packet beyond the linear approximation. To find the first correction 
to the expression for xm obtained before may seem trivial, because if we include the 77-term into 
'ip{K,t), we arrive at an integral that can still be calculated expHcitly. However, as shown in the 
appendix, the correction contains also a contribution of the C-term which is comparable with the 
contribution of the 77-term or greater. After taking into account this contribution we find that 
(see equation (|A-13|) and the comment preceding equation (|A-14p ) 

1 2 

the first correction to a;M ^ —VR'^ > 

K 

where k is a typical interval of k within which lj varies, k = Lu/\dkLu\. If Av is to stay the leading 
correction to vq, the expression on the r.h.s. must be much smaller in the absolute value than 
the correction to xm from which Av has been computed, 

-\vnW'^\A^R\. (21) 

K 

Of course, the condition is relevant only if the effect of higher corrections to xm coming from 
finite a is negligible; that is to say, if the perturbation theory with small parameter a can be 
employed. The corresponding conditions are summarized in equation l|A-lip . It turns out that 
a is small enough if it is much less, at the same time, than l/\/|'7i?J and k. 

For nongaussian packets, write the linear part ofU as 

AKX, 



—00 

where b{K,t) = ae^'^^ and X = x — S^^. The additional correction to Xm, supplementing that 
from which Ai; has been computed, equals the value of X at which the function \U{X,t)\ has 
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maximum. Denote this value by Xm- Since b depends on t, Xm is expected to depend on t, too; 
and since t appears in b through which has the same physical dimension as Xm, the leading 
time-dependent term in Xm is expected to be proportional to ^/ with a coefficient of order 1. 
Thus, if one considers a generic nongaussian packet and wants Av to be the dominant correction 
to V, one must have, in addition to small a, also small |^/|. A term proportional to ^/ is present 
in Xm, for example, if the packet has amplitude (1 + aAK/a)aG with a = const. Such packet 
has two maxima placed symmetrically with respect to the center, located at the points 



Xm = ± 



const H £_j + 0{(t) 

a 



(22) 



On the other hand, it is easy to construct a wide class of packets for which the time-dependent 
part of is exactly zero in the linear approximation. For that purpose, we can use the argument 
of [1^ with exchanged space and time dimensions. Suppose a equals a real nonnegative function 
modulo a phase factor with the phase proportional to AK, a — Ae"?'^^ with A> and q — const. 
Then 

lUunl^ ^ J dKdK'B{K,t)B{K',t) cos[{K ~ K'){X + q)], 

where B{K, t) = Ae^'^^ , and this has obviously maximum at Xm = —q- 

A general analysis of nongaussian packets can be found in the appendix. The estimate of 
the time-dependent part of Xm in equation l|A-12p implies that the conditions under which Aw 
dominates other corrections to v are 

1^/1 and |r;K|a« |A^,j|. (23) 

In addition to that, one must again satisfy conditions l|A-lip in order that the perturbation 
theory with small parameter a is applicable. 

In a general time-dependent medium, the quantities entering the constraints (jA-lip . l(2T|) and 
l(23]) are of order 

. , 1 I L05t 

AU J- ii^-, m-—- 24 

KUJOt K 

From these estimates we find that the perturbation theory is justified only if cr <C {ujSt)~^^^n, and 
Aw cannot be the leading correction to vq for a general packet. However, it becomes the leading 
correction if the packet is Gaussian and has a ^ (aj(5t)~^K. Note that an analogical condition for 
a quantum mechanical particle in an external potential, obtained by flipping over space and time 
directions and identifying wave characteristics with mechanical ones, yields a packet that is much 
larger than the typical distance on which the potential varies. Thus, the post-WKB correction 
to V has no relevance for the description of the motion of a quantum mechanical particle in an 
external potential. 
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For a medium with dispersion relation we must modify the theory sUghtly in order that 
the results of the appendix remain valid. Namely, we must define the phase V' in the geometrical 
optics approximation in a different, and perhaps more natural, way, with omitted logfc. This 
yields an expression for U stripped of the overall factor 1 / \/k, which can be regarded as absorbed 
into a. Using this 'reduced' definition of ip, and assuming V ~ —At^^, we find that the relevant 
quantities are of order 

A^Rr^ k-^At-^, £,1 k-^Ar^, riR^k-'^At-\ (25) 

(The second estimate comes from ipj — — V"/(4fc^).) If we also take into account that k — k, we 
conclude that both sets of constraints l|A-lip and (|23|) are satisfied provided a <^ k, which we 
have assumed from the very start. Thus, for the 'cosmological' dispersion relation, the quantity 
Aw coming from the post-geometrical optics contribution to iji is the dominant correction to v 
for all wave packets whose width in A:-space is small in comparison with the typical value of k. 

3.4 Exact solution 

To support our conclusions, let us construct an exact solution to equation Q. Suppose 

y = -Y2- (26) 

As seen from ^ , V assumes this form for gravitational waves in a dust universe (which expands 
after the law a oc t^). First we find solutions with definite wavelength. If we pass from U to 
U/t and from t to kt, l|10p transforms into the equation for spherical Bessel function with I = 1. 
Moreover, we can see from lfT3| with omitted logarithm that the asymptotics of U for /ct » 1 is 

U = e-«[fc*+i/(fc*)]. This yields 

U ^ -kthf\kt) =(l-^\ e-'^K 



kt^ 

Next we form a wave packet with the amplitude Kac/k. The function U is 

[l + i(p-gs)]e"^''', 

where s — a{x—t),p — l/{kt) and q = a/k. After inserting this into the condition ds\U\'^{sM) — 
we find 

[1-q^ + ip- qsAi f] SM +pq = 0. (27) 

By definition sm = cr{xM — t), therefore the equation for sm determines the worldline of the 
peak of the packet x — XM{t). In fig. [T]this worldline is depicted for fc = 1 and a ~ 1/3. At 
the beginning, the peak of the packet stays within the light cone, but then it moves out. This is 
consistent with the sign Av assumes at late times according to ((HI). 
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Figure 1: Gravitational wave pack;et in a dust universe 
3.5 Spherical packet 

Consider equation ^ with the source 47rj. If we ignore the field V, the solution is 

«(0)(r,i)= J ij(r',i-i?)dV, (28) 

where R = |r — r'|. Suppose the source is nonrelativistic. In the wave zone, u*^*^' can be written 

as 

u("' = i(/ + n.g) + 4n.g, f^fj{r',T)d^x', g = / j(r', r)rW, 

^ j,Z J J 

where t = t ~ r and n = r/r. The magnitude of u^*^-* is 

|u(")| = -(|/|+n.Gi) + ln.G, G = ^Re(/g*), Gi - ^Re(/r). 
r H I/I I/I 

We are interested in the radial coordinate r^^^ at which |m(°) | has maximum for the given direction 
of the propagation of the packet n. To compute this quantity, we expand the derivative 

5,|zi(")| = -- (l/r + n . Gi) - 1 (l/l + n . Gi + n . G 



around tm = t ^ tm, where tm is the value of t at which the function |/| has maximum. In 
fact, we can put r — in the argument of all functions except for |/|', which we must expand 
up to the first order in Ar = r — tm ■ After replacing tm by t we obtain 



Ar' 



(0) 
M 



1 



l/l 



M 



n . GiA^ + i ( \f\M + n . Gi^ + n . G^ 



(29) 



where the index M indicates that the quantity in question is to be taken at r = tm (or, 
equivalently, at r = tm)- By differentiating this with respect to t we finally find that the 
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correction to the group velocity coming from the nonplanarity of the wave packet is 

1 



I/I 



M 



\f\Mt^ 

To estimate the dominant term, note that |/|a4 
packet in the radial coordinate. This yields 



n .GiM+n .Gm) ■ (30) 
' ~l/|A4/o'r, where Gj. is the dispersion of the 



a; 



We can see that the group velocity of the radiation emitted by a compact source is superluminal 
even without a time-dependent medium, merely due to the fact that the factor 1/r shifts the 
maximum of the envelope of the packet the more backwards the closer the packet to the source. 
This is demonstrated in Fig. [H 



SHIFT 
DUE TO 1/r 



CENTER 
OF THE PACKET 



u 




Figure 2: Superluminal velocity of a spherical wave packet 

If we put the field V back into the wave equation, we obtain, in addition to Auf"', another 
contribution to the group velocity Awy. As we shall see, nonplanarity has significant effect on 
the propagation of the packet, but in relation to Auy we can make use of the fact that in a 
sufficiently small domain any spherical wave can be replaced by planar, and identify Awy with 
the contribution to the group velocity Aw computed before. To estimate the relative size of the 
corrections Auy and /\v^°\ consider again, as in section[2l a dust universe. Then Auy ^ ikt)~'^, 
and since the dispersion in x-space ar equals approximately the inverse of the dispersion in 
fc-space (T, we have 

-—f^ - < 1. (31) 

From (pOl we can estimate the next-to-leading term in Aw^^^ as ^v'^^^i ~ WgAw'^''^ where Vg is 
the typical velocity with which the parts of the source move with respect to each other. This 
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yields 

^ ^ (32) 

The leading term in A?j("^ is determined by the form of the packet in the vicinity of the peak 
of the envelope, and can in principle be established by the local measurement of the oscillating 
field. Thus, Awy can be found from local data if it dominates Aw^g^^^, which is the case if the 
velocity of the internal motion of the source is small enough, Vs ^ cr^/fc^. 



4 Conclusion 

When inspecting the equation for gravitational waves in an expanding universe, one immediately 
observes that the group velocity predicted by it is superluminal provided the pressure to density 
ratio is less than 1/3 (which is surely the case in our universe) and the typical wavelength of the 
waves lays well within the horizon. The positive sign of the additional velocity Aw is entirely the 
consequence of the negative sign of the 'external field' V, which in turn refiects the fact that for 
gravitational waves there exists a finite Jeans length of order of horizon length. Note that for 
electromagnetic waves. Aw is exactly zero. The reason is the conformal invariance of Maxwell 
action in four dimensions, which prevents the intensity of the field from having 'tails' inside the 
fight cone; for details, see [3]. To support our calculation of Aw, we have performed a detailed 
analysis of wave propagation in a time-dependent medium. The value of Aw proved correct for 
plane waves, and one would expect that it will be automatically correct for spherical waves, too, 
since group velocity is a local characteristic and spherical waves are locally equivalent to plane 
ones. However, the effect is so small that the effects of nonplanarity turn out to be substantially 
superior to it. For gravitational waves coming from astrophysical sources, the net time shift At 
caused by Aw is at best of order of 10 ps. Such a short time would be surely hard to measure even 
if its beginning and end were defined by single events, which they are not. To determine Ar, one 
should reconstruct the envelope of the pulse in the vicinity of the peak with an error less by many 
orders of magnitude than the period of the pulse; calculate the position of the peak relative to 
the forefront at the moment of the formation of the pulse with the same accuracy; and separate 
out the effects of nonplanarity. A principal obstacle to the first step seems to be that in order to 
carry it out, one needs to know the imaginary part of the oscillating field, which is an auxiliary 
object inaccessible to the measurements. However, the analysis in the appendix suggests that 
the complex field can be reconstructed from the data by a suitable approximate procedure, most 
simply, by determining the Fourier coefficients of the real field seen in the experiment and making 
use of the coefficients with positive /c's only. While this would typically produce a much greater 
error than the measured effect, only the time- independent part of the position of the peak of 
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the envelope (the quantity sq in the notations of the appendix) would be affected. Obviously, 
there is no hope that the outlined procedure to measure At could sometimes be carried out in 
practice. Nevertheless, we believe that the effect is still of some interest since it contributes to 
our understanding of how the propagation of gravitational waves is affected by the expansion of 
the universe. 

Acknowledgement. This work was supported by the grant VEGA 1/3042/06. 



A Propagation of planar wave packets 

A general planar wave packet in a time-dependent medium is described by the integral 

oo 

J a{K)e-'^'^+'^'^''dK, (A-l) 



where 



A^ ^^lj{K,t) -^j = ^AK +-r^AK'^ + -CAK^ + AK = K-k. (A-2) 

2 6 



Introduce the Fourier transform of a, 

CG 

a(s) ^ J a{K)e'^^'dK. 

— oo 

If we write the exponential function in the expression for U as 



where A" = a; — ^, we find 

i 1 
U = a+ -rja" + -Qa'" + . . . , 

where the functions a, a", ... are taken at the (complex) value of the argument s = X . Next, 
we expand the functions a, a", . . . around the point sq at which \a\ has maximum. This yields 

U^Ao+ AiAX + \A2{iri + AX'^) + \a^{C + SiTyAA + AX"^) + ..., (A-3) 

where An is the nth derivative of a at s = sq and AX = X — s^. Suppose a{K) is of the form 
f{AK/a)/a, where |/| has maximum of order 1, f{u) is significantly different from zero only in 
the interval \u\ < 1, and the derivatives of / obey |/*^"-'(u)| < 1 everywhere in the interval |m| < 1. 
For a large class of such functions (for example, for Gaussian functions with the dispersion of 
order a multiplied by appropriately normalized, but otherwise arbitrary polynomials), a{s) is of 
the form /(crs), where f{u) has the same properties as f{u). Then, the coefficients in the sum 
l|A-3P can be estimated as 

A„ = /^")(^so)a"^^", 
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and the sum is in fact a power series in a. The real part of the shift A A" is AX = X — so, where 
X = Xr — X — S^r. Here, as in section [31 we denote by R and / the real and imaginary part of 
the quantity in question. If we pass from AA" to /S.X in the expression for U, we obtain a series 

W = iio + wiAX + u2AX2 + ..., (A-4) 

with each coefficient itself a series, 

1 1 1 

n\ ^ — ' 2 o 

in=0 

(A-5) 

We are interested in the shift of the peak of the wave packet Xm = Xm — (.r- This can be written 
as Xm = So + AXm, where /S.Xm is the value of /S.X at which the function \U\^ has maximum 
at fixed t. Extrema of \U\'^ are zeros of the function 

n 

cx C/o + UiAX + C/aAX^ + . . . , [/„ = ^ (n + 1 - m)(u™<+i_^)ij. 

m=0 

For small enough a we expect \U\ to be close to \a\ aX X ^ sq. Then, the value of AXm can be 
found by solving the equation Ua + UilS.X + J72AX^ + . . . = iteratively, 

2" 



A4.1' = -g-A, Axi5) = A-|A^ a4^;.A-|a^ 



The quantities Uk are generically of order 0-'°+^, except for J7o which is of order a^. Indeed, 

and if we fix, for the sake of simplicity, the overall phase of a so that a{so) is real, we find that Aq 
is real and Ai purely imaginary. From the estimate of Uk it follows that the term proportional 
to A" in the iterative expression for AXm is generically of order cr"~^, so that the procedure is 
meaningful for small enough a. Finally, we expand all terms entering AXm in the powers of a 
to obtain 

AXm = Xmo + Xmi + Xm2 + • • ■ , Xmu c". (A-6) 

Note that sq is generically of order l/cr (the width of the wave packet in x-space) , therefore a 
power series for AXm defines a power series for Xm, too, but starting with a term of order minus 
one rather than zero. If we write 

Uo = A + AA+ Ui^B + AB + ..., C/g = C + . . . , 

where A, B are of order and A^, AB, C are of order a^, we find that the first two terms in 
the expansion of AXm are 

Xmo = ~^ = S, XMi = -^iAA + ABS + C5^). (A-7) 
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(We do not consider the case B = which requires a separate discussion.) The expressions 
for A A, AB and C are obtained from the expansions of first few u's in the powers of ct. A 
straightforward computation yields 

^ = ^0^2/6, B = AoA2r + AIj, AA^-^Ao[A3BiT]i + ^j) + A3n]R] + ^Au>^ 

4A2RVR - A2i{m - el)], AB = [AoAu - AuA2ii)£,i, C = ^(^0^3^ + 3A1/A2/). (A-8) 

Note that only ^ and 77 (the first and second derivative of i/') appear here, although C, (the third 
derivative) is to be included, too, into the formulas for U and \U\^ , if one writes them with the 
accuracy needed for the calculation of AXm up to the term Xmi- The point is that C, enters 
only the term independent on AX in \U\^ and therefore plays no role in the extremization. Note, 
furthermore, that if A21 is zero (the phase of the function a has zero second derivative at the 
maximum of la]), Xmo is zero and AXm is at least of order a. For this class of packets, AXm can 
be made arbitrarily small by shrinking the packet in /c-space. To understand why A21 rather than 
All determines Xmo, note that Au can be adjusted arbitrarily by shifting the packet in fc-space. 
But the shift of the packet by the distance q results only in multiplication oiUun ~ a{X) (the 
function U in the linear approximation) by the factor e**'* , which has no effect on the position 
of the peak of \Uun\ at fixed t. Thus, we can change An without changing Xm in the linear 
approximation; and since the linear approximation is all what is needed to determine Xmo, the 
value of Xmo cannot depend on Au . 



An important special case is the Gaussian wave packet l|T6|). The Fourier transform of a 
Gaussian amplitude is again Gaussian, 



ao = \/27r(Texp ( cr^s^ 



Since this is real, all imaginary A's vanish; and since this is an even function of s, with the 
maximum at the origin, sq = and all odd A's vanish, too. This implies that Xmo — Xmi = 
and the leading term in AXm as well as in Xm is Xm2- For an arbitrary packet with real nonzero 
Ao, A2, Ai, . . . and zero Ai,A^, ... we have 

Uo = {uouDn, = ]:AMC.R + "^^i^r) + ^^26*7/? + ■ • • , C/i = "^iuou^R + juip = ^0^2 + ■ ■ • 
2 

If we insert here the A's for a Gaussian packet, A21, = {—lY{2v — l)!!tT^'', and use the resulting 
Uo and Ui in the formula Xm = —Uq/Ui + . . ., we find 

Xm2= (^^iVr + \<:r^<^^- (A-9) 

Note that the first term can be obtained also from the analytical expression for U in the quadratic 
approximation, with the A/C-expansion of ip truncated after the ?7-term. In this approximation. 
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U can be found from the formula for Uun written down in section [H by replacing 

1 1 

This yields 



where a = 1 — 7//cr^ and A — ^/a? + T^^. The wave packet is Gaussian as in the linear 
approximation, but with a time-dependent width cfx = A/{^/aa) and the peak shifted by Xm = 
^iVR'^^/^- Within the current accuracy, the factor 1/a in the expression for Xm plays no role. 
Thus, Xm effectively equals the first term on the right hand side of (|A-9p as mentioned above. 



Let us find conditions under which our perturbative approach is applicable, and estimate the 
time-dependent part of Xm- We are interested in this part only, because we wish to compare 
different contributions to v — xm rather than to xm- If we use the "reduced" definition of ip, 
with whatever factor in uj depending only on k removed from the term proportional to logw, the 
time-dependent part of Xm is just AXm- (If we left such factor in to, we would have to skip 
from Xm, in addition to sq, also a part of Xmo-) As demonstrated in previous calculations, the 
term Xmu in the expansion of AXm is a sum of products 

mA 

multiplied by coefficients of order 1. For the real and imaginary part of the nth derivative of ip 
we have ^p'^p^^ ~ fm/K"- and t/i^"^ ^ where | — ^k, fj — rjn'^ and k is the scale of k on which 

UJ varies. Thus, if P(^j^^uj) is the sum of terms in Xmu with given vr — ^ i^mR and vi — Yl v^i 
such that n' — 2vr + z// < rt + 1, we have 

P(u.,u.)-i?fl7-^,- (A-10) 

Rewrite this as 



71+1— 



where the first expression is to be used if vr — 0, the second expression is to be used if vi = 0, 
and both expressions are equally good if ^ as well as, vi ^ As can be seen from this 
estimate, we must require 

|r,fl|a2«l, - « 1, (A-11) 

K 

in order that the parameter a can be considered small. This guarantees that all P(ur,v,), as 
well as their sum Xmti, are suppressed the more the higher the value of n. In particular, if all 
three parameters in l|A-lip are of the same order e, P(vn,vi) and Xmu are of order e". In the 
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perturbation theory we can estimate the quantity of interest by the leading term, or terms, of 

the expansion in the small parameter of the theory. For a general packet, the estimate of 

as well as the exphcit expressions for Xmo and Xmi suggest that the leading terms in AXm are 

AXm + Vr<^- (A-12) 

For a Gaussian packet, Xm in the lowest nonzero order in a is 

Xm2 - (^^i + m^''- (A-13) 

Any term in Xm that cannot be obtained from these two by multiplying them by some product 
of powers of the three parameters in (|A-lip . is of the form / n)"^ or ?7_Rcr(7]flO-^)" (is 

proportional either to a product of some powers of ^/ and its derivatives, or to some power of 'qs)- 
We will argue that such terms are actually missing, so that the two terms of the lowest nonzero 
order in a are at the same time the leading terms. The nth term in the iterative expression for 
Xm is proportional to A", that is, to C/q . If we write Uo as 

ran u 

we can see that in order that there exists a term in Xm that assumes one of the two forms 
mentioned above, there must exist at least one expression of the type {'4'even^odd)R which assumes 
the same form with skipped a (since any term in Xm contains at least one such expression). The 
tp^s can be written as 

where 



~ + ?/'^"' V^"'-* + V^"' V^"'^ + ■ • ■ , «^ > 2, 



n, = n. 



The terms in V'n containing only ^7 and its derivatives are of the form i"-^Y^->pj"^^'ip'f'^'' . . ., with no 
factor i contributed by thus the terms of the same kind in ipevenipodd are purely imaginary 
and none of them appears in {ipeveni^odd)R- The quantity ^pn does not include terms containing 
just rjH if n is odd, since there is necessarily at least one odd derivative of ^p present in every 
contribution to in this case, therefore the expression 4'even'4'odd does not include such terms 
either. Thus, there are no terms of the type ^^/k™ or 77^ in {ipeveni'odd)R, q- e. d. As a result, 
the estimate of Xm coincides with the estimate of Xm2- However, one of its two terms can 
be omitted if we note that |^/| < 1/k. (For a general dispersion relation we have ^/ ~ 
and for special dispersion relations we can have, when using the "reduced" definition of ip, also 
1^7 1 < 1/k.) This yields 

Xm - -VR'^^- (A-14) 
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